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ABSTRACT 

We present a study of three-dimensional structure of the Small Magellanic Cloud 
(SMC). The V- and /-band light curves of the fundamental mode RR Lyrae stars 
(RRab) obtained by the Optical Gravitational Tensing Experiment (OGLE)-III 
project we re utilized in order t o com prehend the SMC structure. The lFe/JI]—P—(/) 3 i 
relation of lJurcsik fc Kovacsl l|1996ll is exploited to obtain the metallicities. From the 
three-dimensional RRab distance distributions, northeast (NE) arm and main body 
of the galaxy is identified. Combining metallicities with spatial distribution of these 
tracers, no radial metallicity gradient in the SMC has been detected. Dividing the 
entire sample into three parts: northeastern (NE), central and southwestern (SW), 
we hnd that the central part has a significantly larger line of sight depth as com¬ 
pared to rest of the parts, indicating that the SMC may have a bulge. Results ob¬ 
tained f rom the /-ban d data seem to be reliable and were further substantiated us¬ 
ing the ISmoled (120051) relation. Distribution of SMC RRab stars were modeled as 
a tri-axial ellipsoid. Errors in structural parameters of the SMC ellipsoid were ob¬ 
tained from Monte Carlo simulations. We estimated the axes ratios of the galaxy 
as 1.00 ± 0.000 : 1.310 ± 0.029 : 8.269 ± 0.934, the inclination of the longest axis 
with line of sight i = 2°.265 ± 0°.784, and the position angle of the line of nodes 
^lon = 74°.307 ± 0°.509 from the variance weighted /-band determinations. 

Key words: stars: variables: RR Lyrae-stars:fundamental parameters - stars: Popu¬ 
lation II - galaxies: statistics - galaxies:structure - galaxies:Magellanic Clouds 


1 INTRODUCTION 

The Small Magellanic Cloud (SMC) is one of our nearby 
satellite galaxies. It is a gas-rich, dwarf irregular galaxy 
connected by a hydrogen gas and stellar bridge to the 
Large Magellanic Cloud (LMC) and located at a distanc e 
of around 60 kpc (IWesteriundlll997l : ICraczvk et al.ll2014h . 
Like LM C, the SMC also exhibits a bar which is less pro¬ 
noun ced llWesterlundl Il997l : ISubramanian fc SubramaniamI 
I 2 OI 2 I ). The complexity in the structure, dynamics and evo¬ 
lution of the SMC might be attributed to the influence of 

J ravitational interac tions with the LMC and the Galaxy 
Putman et al.lfl9^ ). The SMC has a relatively lower heavy 
metal abundance, lower dust content and significantly larger 
line of sight depth as compared to the LMC and the 
Galaxy (IStanimirovic et al.||l99 ^: IStanimirovic et al.l l2004l : 

ISubraimanian fc SubramaniamI I 2 OI 2 I) . With the prolifera- 
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tion of accurate and precise variable star data generated 
from modern astronomical surveys and space missions along 
with their well-established theoretical understandings have 
helped us to gain better knowledge about distances of as¬ 
tronomical objects and resolved many issues related to dis¬ 
tances and distance related parameters such as Galactic 
and e x tragalactic structures (ISubramanian fc Subramaniam 


2012; |Ka j3akos_&^atzidimitriorJ l2012l : IPietrzvski et al 


2012; iGraczvk et al. 2014 ). For instance, to gain a bet- 


ter knowledge of t he geometric structure of the SMG, 
IGraczvk et al.l (|2014| ) determined the distance to the SMG as 
(m — M) = 18.965±0.025(stat.)±0.048(sys.) mag which cor¬ 
responds to a distance of 62.1 ± 1.9 kpc using detached, long 
period eclipsing binaries. Well-detached eclipsing binaries 
serve as a source of accurate and precise distance indicators 
and help in the calibration of the zero point of the cosmic dis- 
tance scale with an accuracy of about 2% dPietrzvski et al.l 

I 2 OI 3 I) . 


The long-standing and highly controversial questions 
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of the SMC’s three-dimensional structure and depth along 
the line of sight have been greatly facilitated through 
the identification and characterization of a statistically 
significant number of different stellar tracers in vari¬ 
ous large photometric surveys. The uniform data sets 
obtained for a large number of different stellar tracers 
obtained in these surveys in the SMC provide a unique 
opportunity to study the characteristics of the host 
galax y to a g reat detail ilSu bramanian fc S ubramaniamI 


2OI2I: 


Kapakos fc Hatzidimitrioul 


Subramanian fc SubramaniamI 2014) . 


2OI2I: 


The accurate determination of the geometrical param¬ 
eters of the SMC plays an important role in modeling the 
interactions between the galaxy, LMC and the SMC. In this 
study, the ‘standard candles’ such as the RRab stars have 
been chosen to study the various parameters of the SMC. 
RRab stars are the RR Lyrae stars which pulsate in the fun¬ 
damental mode. RR Lyrae stars are metal-poor, low-mass 
and core-h elium burnin g stars that undergo periodic radial 
pulsations (ISmithll2004 . They serve as excellent tracers of 
the oldest population of stars in a host galaxy. Their con- 
stant mean lumin osity makes them ideal standard candles 
dPritzl et al.ll2oiTh . They serve as invaluable resource which 
help in exploring the metallicity and structure and hence 
provide substantial insight into the un derstanding of the 
formation and evo lution of the galaxy (IParisi et al.l I2OO9I : 
iDeb fc Singh|[20l4 . They have been used as proxies of the 
true, three-dimensional spatial distribution of host galaxies. 
Their presence in large numbers in different galaxies has led 
to an improved understanding of various geometric param¬ 
eters of these galaxies. Well-sampled light curves covering 
the entire pulsational period of the RRL stars allow deter¬ 
mination of their properties with high precision. The derived 
structural parameters extracted from their light curves can 
be used to determine metallicity, interstell ar extinction, dis¬ 
tance and spatial distribution of the stars (IPietrzvnski et ahl 

l2012h . 

are a few studies in the literature 
determining the geometrical struc- 
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Kapakos fc Hatzidimitrioul l2012l: iHaschke et al.l 


199£: 


2012 : 


2012 ; 


Subramanian fc Subramaniamll2014l~). In v estiga ting 12 pop¬ 


ulous clusters in the SMC, Crowl et al.1 (l200lll determined 
the axes ratios approximately as 1 : 2 : 4 modeling the SMC 
as a triaxial galaxy with declination, right ascension and 
line of sight depth as the three axes. Using the red clump 
and RR Ly rae stars in the V and J-band d ata ta ken from 
OGLE-III, ISubramanian fc SubramaniamI ll2012ll derived 
various geometrical parameters of the SMC. Their study 
was based on the relative position of different regions of the 
SMC with respect to the SMC mean distance. The param¬ 
eter determinations of the SMC in their work do not quote 
the errors. It is to be noted that the geometry of the SMC 
is a function of the de-projected Cartesian coordinates as 
well as their errors. Accurate parameter values can only be 
determined when the errors are incorporated in projected 
Cartesian coordinates along with the reliable distance 
determinations. Also, the metallicity determination of the 
RR Lyrae stars were not attempted in their work which is 
very crucial in understanding of substructure of the galaxy. 
The axes ratios, inclination angle of the longest axes along 


the line of sight (i) and position angle of the line o f nodes 
(dion) obtained by ISubramanian fc SubramaniamI (l2012|j 
from the distribution of the SMC RR Lyraes were found 
to be 1 : 1.33 : 6.47, 0°.4 and 74'’.4, respectively from the 
entire dataset. Taking different coverage of the dataset and 
removing the isolated northwestern fields, different struc¬ 
tural parameters of the SMC ellipsoid were also at t empte d 
and estimated by ISubramanian fc SubramaniamI ll2012h . 
On the othe r hand, using the 454 OGLE-III RRab data in 
the U-band. lKapakos fc Hatzidimitri^ (l2012l l used the line 
of sight distances and the metallicity values to investigate 
the possible presence of different structures containing 
different populations in the SMG. Using spherical shells 
of different radii of 2.5, 3.0 and 3.5 kpc, they obtained 
the axes ratios of the SMG ellipsoid as: 1 : 1.21 : 1.57, 
1 : 1.18 : 1.53 and 1 : 1.23 : 1.80, respectively, whereas, 
inclination angle and position angle of line of nodes of the 
SMC were not calculated in their analy ses. From an analysis 
of radial velocities of 150 carbon stars, iKunkel et al.l (l2000lj 
estimated the inclination of the SMC orbital plane to 
be i = 73 ° ± 4° relative to t h e plan e of t he sky. Using 
Cephe ids, ICaldwell fc CoulsonI lll986h and iGroenewegenI 
ll2C)OOll estimated i as 70° ± 3° and 68° ± 2°, respectively. 
The positi on angle of line of nodes was estimated to be 
~ 148° dGroenewegenl I 2 OOOI L Estimates of inclination 

angle and po si tion angle of line of nodes obtained by 
iHaschke et :ff] ll2012ll from the I-band OGLE III SMC 
RRab data are 7° ± 15° and 83° ± 21°, respectively. 
From the study of distance distrib ution of Ceph e ids in 
the SMC using the OGLE-III data, IHaschke et al.1 (l2012|j 
found that these young populations are differently oriented 
than the old population RR Lyrae stars. The inclination 
angle and position angle parameters of the SMC obtained 
from Cepheids by IHaschke et al.l 1 I 2 OI 2 II are 74° ± 9° and 
66 ° ± 15°, respectively. Other recent estimates of inclination 
angle and position angle of lines of nodes for the SMC 
obtained froi n Cepheids are 64°.4 ± 0°.7 and 155°.3 ± 6°.3, 
respectively (ISubramanian fc Subramaniamll2014^ . 

In this paper, we use available V and I-band OGLE- 
III data in order to determine the absolute magnitudes and 
metallicities of RRab stars. These measurements along with 
equatorial coordinates, viz., right ascension (a) and decli¬ 
nation ((5) have been utilized to get an insight into the un¬ 
derstanding of the three-dimensional structure of the SMC 
and its metallicity distribution. The availability of data in 
two different bands and the use of different methodologies 
of distance determinations provide a unique opportunity to 
compare and contrast various structural parameter determi¬ 
nations. The roadmap of the present investigation is to use 
these tracers for the distance determinations and obtain dis¬ 
tances relative to the center of the SM C. The principal axis 
transformation method as described in iDeb fc SinghI (l2014fl 
was used to get the viewing angle and geometrical parame¬ 
ters of the galaxy. Section [5] describes the data used for the 
analysis in the present study. Fourier decomposition method 
and sample selection criteria are described in section |3] Sec¬ 
tion 2] focusses on extinction measurements for the OGLE- 
III data a long the l ine o f sight using the SMG extinction 


mao from IZaritskvl (Il999l'); Zaritskv et al. 

(|200dL Determi- 

nation of metallicities llJurcsik & Kovact 

1996(1 and abso- 

lute magnitudes (iKovacs & Walked 200ll; 

Catelan & GortesI 


I 2 OO 8 II from U-band and /-band data calibrated to the V- 
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band are described in section [5] Viewing angle parameters 
and geometrical parameters of the SMC determined for the 
above dataset are described in section [6l Line of sight depth 
and existence of metallicity gradient are examined in sec¬ 
tions 0 and [51 respectively. Results obtained from the above 
anal yses were furthe r substantiated as described in section[9] 
with ISmoled (l2005l l metallicity relation using the /-band 
light curve data and exp loiting the Mv — [Fe/H] relat ion of 
I Caftan fc Corti3 (|2008l L Exponential disk model an d lKind 
(Il962tl three parameter model for the radial number density 
distribution of RRab stars in the /-band are discussed in 
section fTOl Lastly, summary and discussions of the present 
investigation are laid down in section [TTI 


2 THE DATA 

RRab stars for the present study were selected from OGLE- 
III catalog of variable stars that consists of 8-year archival 
data identified and characterized by the Fourier coefficients 
of the light curves llSoszvhski et al.ll2009ll . The catalog com¬ 
prises 1933 RRab stars in the /-band having a mean pe¬ 
riod of < Pab >= 0.596 days wit h typical accuracy of 
apjP ~ 10~® llSoszviiski et alJl2009l L The catalog is based 
on the observations carried out at the Las campanas Ob¬ 
servatory, Chile with a 1.3-m telescope. The time series ob¬ 
servations of the SMC consist of 800 nights between June 
2001 and May 2009. The OGLE field in the SMG covers 
nearly 16 deg^ covering the bar and the wing. Most of the 
observations were carried out using the Gousins /-band filter 
with exposure time of 180 s having an average of 400 photo¬ 
metric observations. The catalog also contains V-band light 
curves of 1887 stars having an average of 30 data points per 
light curve with the exposure time of 225 s (ISoszvhski et al.l 
l2009lL 


3 FOURIER DECOMPOSITION METHOD 
AND SAMPLE SELECTION 

The Fourier decomposition method was used to obtain V 
and /-band light curve parameters of the SMC RRab stars. 
The li ght curves were fitted with a Fo urier sine series of the 
form (iDeb fc Singhll2009l . l2010l. |2014^ 

N 

m(t) = Ao+ '^Ai sin [iuj{t - to) + (t)i] , (1) 

i=l 

where m{t) is the observed magnitude, Ao is the mean mag¬ 
nitude, a;=27r/P is the angular frequency, P is the pulsa- 
tional period of the star in days and t is the time of obser¬ 
vation. to represents the epoch of maximum light (either in 
the V or the /-band) and is used to obtain a phased light 
curve which has maximum light at phase zero. Ai’s and ifii’s 
are the ith order Fourier coefficients and N is the order of 
the fit. Eqn.[T]has 2N + 1 unknown parameters. To solve for 
these parameters, at least the same number of data points 
are required. The light curve phase was obtained using 

P \ P 

Here 4? € [0,1] represents one pulsation cycle of the RRab 
stars. The pulsation periods (P) and the epoch (to) are taken 


from the O GLE catalog. The F ourier parameter code as 
described in jPeb fc SinghI ll20ldl was used to obtain vari¬ 
ous Fourier parameters on the right hand side of Eqn. [T] 
Fifth order and seventh order Fourier fits were employed to 
model the RRab V and /-band light curves, respectively. 
The Fourier fitted V and /-band light curves of a sam¬ 
ple of RRab stars are shown in Fig. [T] The phase differ¬ 
ences, (f>n = 4>i — and amplitude ratios, Rn = {Ai/Ai), 
i > 1 were evaluated and standard errors were determined 
using the formulae of iPeb fc SinghI (|2010|L A clean sample 
of RRab stars for the present analysis were selected with 
the criteria based on OGLE-determined periods, the mean 
magnitude (Aq) and the peak-to-peak V and /-band am¬ 
plitude {Av,Ai) determined from the Fourier analysis of 
the phased light curves. RRab stars with periods P ^ 0.4 
days, mean magnitude Ao 18.0 mag (/-band)and ampli¬ 
tude 0.1 ^ A/ ^ 1.2 mag were chosen for the /-band light 
curves. On the other hand, to select the U-band light curves, 
the same criteria on periods has been applied. While the fol¬ 
lowing criteria were applied on Ao and amplitude: Ao ^ 18.5 
mag (U-band) and amplitude 0.2 ^ Ay ^ 1.5. These selec¬ 
tion criteria were applied in order to collect a clean sample 
of RRab stars of the SMC which are free from any contami¬ 
nation due to the foreground objects of our Galaxy and the 
LMC. The application of these selection criteria along with 
furthe r application of compatibility test of jjurcsik fc Kovacsl 
(Il996ll reduces the number of RRab stars in the V and I- 
band light curves to 1142 and 1543, respectively for the phys¬ 
ical parameter estimation and dete rmination of struc t ure o f 
the SMC. The compatibility test of ijurcsik fc Kovacsl (Il996fl 
ensures to reject those RRab stars that have deviation pa¬ 
rameter values Dp > 5. 


4 INTERSTELLAR EXTINCTION 

Interstellar medium plays an important role in dimming 
and reddening of the star light. In order to determine ac¬ 
curate distances to stars, interstellar reddening must be ac¬ 
counted for. Various methods to determine interstellar ex¬ 
tinction exist in the literature using different tracers. Ob¬ 
servations of multiband photometry available for RR Lyrae 
stars se rve to resolve internal extinction in the case of host 
galaxy (jPeicha fc Stane^l2009l l . For the determination of in¬ 
terstellar extinction in the case of individua l RRab stars in 
this p aper, we adopted the IZaritskvl lll999ll : IZaritskv et al.l 
(|2002| . hereafter, Zaritsky map) extinction map for the SMC 
which comes in FITS (Flexible Image Transport System) 
format. The SMC extinction map for cool stars is down¬ 
loaded from Q, since the cool star data are characteristics of 
the RR Lyrae with Tefi ~ 6500 K. The Magellanic Cloud 
Photometric Survey (MCPS) contains UBVI photometry 
and e xtinction map for t he central 18 deg^ area of the 
SMC (IZaritskv et al.l[2002l L For a given RRab location, the 
map returns the U-band extinction value. Ay- The FITS 
image of the Zaritsky map containing the extinction val¬ 
ues (A„) is an array of (240 x 280) pixels with reference 
pixel coordinate as (X, Y) — (120,140) corresponding to the 
{RA,DEC) = (12°.750, -72°.699). The RA and DEC pixel 
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Figure 1. A collection of randomly selected sample of RRab light curves with their corresponding OGLE IDs and deviation parameters 
(Dp) are shown. Left panel shows the phased light curve data in R-band while the right panel depicts the same in the 7-band. 
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Figure 2. Two-dimensional color bar plot of the interstellar extinction values {Ay) and extinction corrected mean magnitudes (Vq) of 
RRab stars. 
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scales in the image are ~ 0.0167°/pixel which correspond to 
I'/pixel. The IRAF task ‘imhead’ was used to get the de¬ 
tails of the FITS image header. The number of stars used 
to derive extinction for map pixels wer e at least 3 whe re the 
median extinction value was adopted llZaritskvlll999l l. The 
Zaritsky map has regions of high and low spatial resolution 
but relatively uniform signal-to-noise ratio. The details of 
smoothing algorithm applied to d erive the extinct ion values 
for the map pixels are outlined in IZarits"^ lll999ll . The fol¬ 
lowing procedure was adopted to get the extinction values of 
the sample of RRab stars. RA(q:) and Dec(5) of the sample 
stars were first converted into X/Y pixel coordinates using 
sky2xy followed by running getpix on the FITS image to 
get interstellar extinction {Av) at that X/Y position. Both 
sky2xy and getpix are the codes from WCSTOOLS 0. For 
a few stars, we get a null value of Av- For these stars, the 
mean value of Av = 0.155 mag appropriate for the SMC was 
used. The average value of Av for the OGLE-III RRab’s is 
0.155±0.127 mag which translates into a mean E{B — V) of 
0.048±0.039 mag. This is comparable to the mean E{B — V) 
value of 0.054 r nag and E(V — I) va lue of 0.07 ± 0.06 mag, 
respectively by ICaldwell &: Cou pon (Il 985ll obtained from 
the analysis of 48 Cepheids and Haschke et al. from 

the analysis of 1529 SMC RRab stars. Since no formal er¬ 
rors are given by the map for the individual stars, an error 
margin of 20% is assumed in the individual extinction val¬ 
ues if these estimates from the map are taken to be reliable. 
The extinction maps derived from the OGLE RRab stars by 
iHaschke et al.l ll201ll . hereafter, Haschke map) were found to 
be in good agreement with those obtained from the Zarit¬ 
sky reddening map. Fig.[2]depicts two-dimensional color bar 
plot of interstelar extinction values [Av) and extinction cor¬ 
rected mean magnitudes (Vb) of RRab stars, x and y in the 
figure denote projected Cartesian coordinates. 


5 METALLICITIES AND ABSOLUTE 
MAGNITUDES 

Metallicities of a large number of stars found in a certain 
galaxy play pivotal role in understanding the formation of 
various structures of that galaxy. Assembly of stars hav¬ 
ing certain metallicity ranges concentrated in a particular 
part of a galaxy tells the epoch of its formation. Although 
the metallicity is an important parameter for understand¬ 
ing chemical properties of the underlying stars and hence 
that of the distant host galaxy and its formation, it is hard 
and time-consuming to obtain the metallicity spectroscopi¬ 
cally for a large number of stars. Calibrated relationships 
between various photometric light curve parameters and 
spectroscopic metallicities of RRab stars provide efficient 
and robust ways for determining the metallicities in var¬ 
ious globular clusters, the Galaxy, the SMC and the LMC 
by recent ground-based automated surveys. Previous studies 
have shown that metallicities of RRab stars can be derived 
from an empirical relation connecting spectroscopic metal¬ 
licitie s and the light curve parameters as p r oxies llSimon 
W8S; lJurcsik fc Kovacj 1 19961 : ISandaeel l2004l : iNemec et al. 
201311 . These studies have shown that metalicities of RRab 


^ http://tdc-www.harvard.edu/wcstools/ 


stars can be linked the period with some intermediate pa¬ 
rameters such as amplitude (A), Fourier phase parame¬ 
ters < 621 , /> 3 i lJurcsik_^Kov^ 1 19961 : lAlcock et al.l 120001 : 
ISandagjl20()4l : iNemec et al.l 2013 1. The correlation between 
the light curve structures of RR Lyrae stars such as the 
period (P), Fourier coefficients and their metal- 

lici ties wa s first noticed and studied in the pioneering work 
by[s imoni lll988ll. This work h a s mot ivated many investiga- 
tors such as Jurcsik fc Kovac'3 lll996ll and his collaborators, 
ISandaeel (l2004ll to derive semi-empirical relations involv¬ 
ing light curve structures and spectroscopically determined 
metallicities. 

The unprecedented amount of photometric light curve 
data generated from the recent automated surveys such as 
OGLE, AS AS (All Sky Automated Survey), and highly ac¬ 
curate and precise stellar data obtained from the NASA’s 
Kepler mission have provided a new avenue into the under¬ 
standing of various astrophysical problem s more accurately. 
Altho ugh the [Ee/H] —P — (j> 3 i relation of ijurcsik fc Kovacsl 
lll996ll is widely used in the determination of metalicities 
of RR Lyrae stars from their photometric light curves, ac¬ 
curate light curve data of these stars generated from the 
NASA’s Kepler mission have helped in refining the model 
parameters of existing \Fe/H] — P — (/31 rela tion and cal¬ 
ibrat e new relations based on the new data llNemec et al.l 
l2013ll . These relations will prove to improve our knowledge 
in understanding the chemical history of the stellar tracers 
as well as structure of the host galaxy from their distance 
distributions. Metallicities [Pe/JLI spec of the corresponding 
light curves in iNemec et al.l (120131) were obtained from the 
high resolu tion spectroscopic meas urements.The empirical 
relation of ljurcsik fc Kovacsl lll996ll connecting (/si in the 
V- band, period P and the metallicity [Fe/H] is given by 

[Pe/Pljif =-5.038 - 5.394 P-t 1.345 (/ 31 , a = 0.14. (2) 


On the other hand, the non-linear relation connecting 
P — </ 3 i — [Pe/P] using 26 RRab stars obtained by the 
NASA’s Kepler sp ace telescope in the Ap-band is given by 
(INemec et al.|[2013 ~) 

[Pe/P]Ni3 = (-8.65 ± 4.64) -b (-40.12 ± 5.18)P -b 
(5.97 ± 1.72)P(/3i -b (6.27 ±0.96)((/3l)^ a = 0.084 (3) 


where (/li is the mean 4>3i{Kp) value and [Fe/H] is the 
spectroscopic metallicity value obtained from t he recent new 
high resolution spectroscopic measurements dNemec et al.l 
I2OI3I). In order to apply the N13 relation to the U-band 
data, INemec et al.l (I2OIIII gave the following inter-relation 

CaiiV) = CsiiKp) - (0.151 ± 0.026). (4) 

The above inter-relation was derived bv iNemec et al.l (l2013fl 
using the AS AS U-band and the Kepler photometric light 
curve based on only 3 stars. The inter-relation was recently 
refined using V band photometric data of 34 common stars 
in the Kepler fiel d and the following relation was obtained 
djeon et al.ll2014h 

= </3i(Pp) - (0.174 ± 0.085). (5) 


In t he present analysis , we h ave determined metallicities us- 
ing JuTcsik__fc_Jioyacj dl996l l linear relation. Although the 
INemec et al.l d2013l l relation is based on the highly accu¬ 
rate and precise Kepler data, we have found from numerical 
analysis that formal errors on the calculated [Fe/H] values 
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Figure 3. Left panel shows LMC and SMC RRab stars in galactic coordinates (1, b). Right panel depicts polar plot of distance distribution 
of LMC and SMC RRab stars in galactic coordinates {l,b). The dashed lines are set at distances D = 25, 50, 60, 70 kpc, re spectively. 
LMC and SMC stars are shown as red and blue data points, respectively. LMC RRab stars are taken from iDeb &: SinghI 1120141 . The tails 
of distance distributions of the two Magellanic Clouds account for possible contaminations from the foreground and background RRab 
stars. 


from the derived relation turn out to be very large using 
the pr opagation of error formula of iBevington fc RobinsorJ 
(l2003l l. This might be attributed to the larger errors in cj> 3 i 
values for the OGLE-III data and more number of complex¬ 
ity parameters of the derived relation. Hence because of lim¬ 
itations of the relation restricted to highly accu r ate da taset, 
we have to leave out using the iNemec et al.l (l2013l i rela¬ 
tion for the OGLE-III dataset. The absol ute magnitudes are 
deter mined using the following relation (iGatelan fc Cort^ 

l2008h 

Mv = (0.23 ± 0M)[Fe/H] + (0.984 ± 0.127), (6) 

where [Fe/H] is the metallicity in the lZinn fc Wes3 (Il984l i 
scale. The metall i city v alues obtained from Eqn. [2] are in 
ljurcsik fc Kova^ lll996li sc ale, which can be t ransformed 
into the met allicity scale o f IZinn fc Wesj (Il984l l using the 
relation from ljurcsiki (Il995l i: 

[Fe/H] = [F'e/H ]jk -0.88 ^ ^ 

^ 1.431 ^ ’ 

Determination of absolute magnitudes of the present 
sample of RRab stars using Eqn. [6] is denoted as ‘MVC’. 
The absolute magnitudes of the selected samples of RRab 
stars were also determined using the following empirical re¬ 
lation in term s of the period (P) and the Fourier coefficients 
{A\ and As) (iKovacs fc WalkCTl200ll : I Arellano Ferro et al.l 

I 2 OIOII : 

Mv = -1.876 log P- 1.158Al -k O. 82 IA 3 + 0.41. (8) 

Absolute magnitude determination using this relation is re¬ 
ferred to as ‘MVF’. Both the methods of absolute magni¬ 
tude determination were applied in V and 7-band to cal¬ 
culate distances of the present sample of RRab stars. Suit¬ 
able inter-relations between the Fourier parameters in the 


V and J-ba nd were utilized as discussed in IPeb &: SinghI 
Geometrical parameter determinations of the 
SMG using two different relations in two different bands car¬ 
ried out in this paper will help us in improving their esti¬ 
mates. Fig.[2]shows polar plot of distance distribution of the 
SMG RRab stars (7-band calibrated to E-b and) in galac¬ 
tic (L 6) coordinates. LMG RRab stars as in IPeb k, SinghI 
11201411 are also overplotted. The two-dimensional density 
cont ours of the dis t ribut ion of the LMG RRab stars taken 
from IPeb fc SinghI (l2014l 'l and the SMG RRab stars in the 
present study are shown in Fig. 2) The star symbol denotes 
the location of the centroid of the present sample. Two- 
dimensional color bar plot of distance distributions of SMG 
RRab stars is shown in Fig. [3 


6 STRUCTURE OF THE SMC 

Availability of a statistically significant sample of RRab 
stars from the OGLE-III data in 7 and E-band provides 
an opportunity to compare and contrast the two different 
bands to geometric parameter determination of the SMG. 
Furthermore, structural parameter determinations of the 
SMG substantiated using two independent distance deter¬ 
mination formula will also help us in improving the struc¬ 
tural parameters of the SMG. In order to determine the 
structure, we have converted the right ascension (a), dec¬ 
lination (d) and the distance {D) into the corresponding 
Cartesian coordinates {x,y,z). Let us consider the Garte- 
sian coordinate system {x, y, z) which has the origin in the 
center of the SMG at {a,S,D) = (ag, So, Do). The z axis is 
pointed towards the observer, x-axis is antiparallel to the 
a-axis, the y-axis is parallel to the (5-axis. Do is the dis¬ 
tance between the center of the SMG and the observer. 
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Figure 4. Left panel shows two-dimensional density contours of the LMC RRab stars taken from [Peb &: SinghI l[20l3). Right panel 
shows the same for the SMC RRab stars in the present study. The star symbol denotes the location of the centroid. 
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Figure 5. Two-dimensional color bar plot of distance distribution 
of SMC RRab stars as a funtion of {x, y) coordinates. 


D is the observer-source distance. (ao,^o) are the equito- 
rial coordinates of the center of the SMC. The center of 
the SMC is taken as (ap.^o) = (0^5 3”"31°. -72°59'15".7) 
dSubramanian fc SubramaniamI l2012fl . The {x,y,z) coor- 
dinates are obtained using the transformation equatio ns 
llvan der Marel fc Cionill200ll : IWeinberg fc NikolaevfboO J l: 

X = —D sin(a — ao) cos 5, 
y — DsinS cos So — D sin So cos (a — ao) cos S, 
z = Do — D sin <5 sin So — D cos 5o cos a — ao cos S. 

The coordinate system of the SMC disk {x',y',z') is the 
same as the orthogonal system {x,y,z), except that it is 
rotated around the 2 -axis by the position angle 8 coun¬ 


terclockwise and around the new a;-axis by the inclina- 
tion angle i clockwise. The coordinate transformations a re 
llvan der Marel fc Cioni|[200ll : I Weinberg fc Nikolaevllioolh : 


x' 


cos 8 sin 8 0 


X 

y' 

= 

— sin 8 cos i cos 8 cos i — sin i 


y 

z' 


— sin 8 sin i cos 8 sin i cos i 


z 


The errors {i7x,o-y,ax) in the Cartesian coordinates were 
obtained following th e pro pagation of error formula 
llBevington fc RobinsonI l2003li . Observed distribution of 
RRab stars in the SMC was modeled by a triaxial ellipsoid. 
Properties of the ellipsoid are obtained from the moment 
of inerti a tensor using the principal axes transformation 
method (IPeb fc Singhll2014h . The axes ratios where i = 
0,1, 2, inclination of the longest axis along the line of sight 

(i), position angle of line of nodes (0ion) along with their as¬ 
sociated errors were calculated using the Monte Carlo simu¬ 
lations. The Monte Carlo simulations for finding errors in the 
geometrical parameters were obtained with the assumption 
that the probability distribution of errors in the parameters 
are Gaussian and the errors in the obtained (x, y, z) values 
are reliable. The observed projected Cartesian coordinates 
[x, y, z) were simulated with the following steps: 

(1) The observed Cartesian coordinates {ri,ari) were ob¬ 

tained from the transformations connecting them with the 
given {ai,Si,Di) alon g with the application of p ropaga- 
tion of error formulae feevington fc Robinsonl[2003h . where 
n = {xi,yi,Zi) and i = 1,2,... ,N and 

N is the number of data points. 

(2) Vi’s were taken as the centroid of the Gaussian. 

(3) Box-Muller method was then used to get Gaus¬ 
sian distributed points having ar^s as errors and Vi’s as 
centroids. Essentially, this c onsists of the following steps 
llBevington fc Robinsonll2003l l: 

ri =ri + Gnki, 

with ki = 2.0 log (ui) cos (2.07rt;i), 
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where {ui,Vi) are pairs of uniformly distributed random 
numbers. 

Principal ax is tran sformation method as described in 
iDeb fc SinghI (I20l4) was then applied to the simulated 
Cartesian coordinates of the SMC RRab stars to obtain 
viewing angles and other geometric parameters of the galaxy. 
The Monte Carlo simulation was carried for 10® steps and 
corresponding values of the parameters were determined in 
each iteration. The distribution of various geometric param¬ 
eters were then obtained after binning with a proper binsize. 
Distribution of the geometric parameters were found to be 
Gaussian. Hence, three parameter Gaussian fitting was ap¬ 
plied to the histogram distributions. The o values of the 
fitted Gaussian distributions were taken as errors in various 
parameters. Viewing angle parameters (i, 0ion) and geomet¬ 
rical parame ters such as axes ra tios (S'o/S'o, ^i/So, S 2 /S 0 ) 
described in iDeb k, SinghI ||2014| ~) obtained from the entire 
dataset of SMC RRab stars are shown in Figs.|6]and[3 Their 
values are listed in Table [T] 

From the Cartesian three-dimensional distribution of 
RRab stars as shown in Fig.| 8 l it is quite evident that north 
east (NE) arm is nearer to us than the SMC main body. 
In order to find exact location of NE arm of the SMC, we 
have used an B n statistical package ch a nqepoint for change - 
point analysis (IKillick fc Eckle'^l2014l : iKillick et al.l[20l^ . 
The function cpt.meanvar() was used to find changes in 
the mean and variance for z distribution. Locations of main 
body and NE arm of the SMC were thus identified from the 
changepoint analysis. Principal axis transformation method 
was then applied to the SMC main body since its distribu¬ 
tion resembles a three-dimensional ellipsoid. This is done in 
four cases (Two in each bands corresponding to the methods 
‘MVF’ and ‘MVC’). Results obtained from the principal axis 
transformation method are shown in Table [2] On the other 
hand, NE arm of the SMC resembles like a plane. Therefore, 
plane-fitting procedure was applied to find its viewing angle 
parameters such as inclination and position angle of line of 
nodes. NE arm of th e SMC was fitted wit h a plane-fitting 
equation of the form (iNikolaev et al.|[2004h 

Zi = c-\- axi -\-hyi, i = (9) 

where N is the number of data points. The inclination angle 
(i[°]) can be obtained from 

( 1 ^ 

and if we define 7 = arctan(|a|/|t(|), then 


7 


if a 

< 

0 

and 

b 

> 

0 , 

7 


if a 

> 

0 

and 

b 

< 

0 , 

7T 

2 

+ 7 

if a 

< 

0 

and 

b 

< 

0 , 

7T 

2 

+ 7 

if a 

> 

0 

and 

b 

> 

0 , 

0 


if a 

= 

0 

and 

b 

/ 0 , 

sign(a)f 

if 

a / 

0 and 

b 

= 0 , 

undef. 

if 

a 


5 = 0 . 





A weighted plane-fitting procedur e using mpfitfun in ID L 
was used to fit NE arm of the SMC llMarkwardtll200^ . l2012ll . 

http://www.r-proj ect.org/ 





Figure 7. Distribution of S 2 /S 0 ratio for 10® Monte Carlo iter¬ 
ations in histograms. Solid red lines denote Gaussian fit to the 
resulting distributions. 


Results of the plane-fitting procedure have been listed in 
Tabled Errors in i and 9ion were obtained from the prop¬ 
agation of error formula. From Table [2] one can find that 
the values of i are comparable to each other within the er¬ 
ror bars indicating that the values obtained using two dif¬ 
ferent datasets and empirical relations are reliable. On the 
other hand, if we compare the values of viewing angle pa¬ 
rameters, then we find that Sion values of the SMC main 
body are significantly different for the two bands. Difference 
in the Sion arises due to poor sampling in distance distri¬ 
bution in the case of V-band as compared to those in the 
7-band. We take the results obtained from the 7-band data 
to be more reliable because of larger number of stars and 
better phase coverage of light curves as compared to the 
corresponding V-band light curves. The viewing angle pa¬ 
rameters obtained for NE arm of the SMC indicate that 
it is slightly misaligned with the SMC main body. Cumu¬ 
lative distribution of extinction corrected magnitudes (Vo) 
and metallicities ([7^e/77]) of NE arm (sample 1) and SMC 
main body (sample 2) are shown in Fig. In order to test 
whether the two samples came from the same distribution at 
a significance level of 0.0 5, we use two samp le Kolomogorov- 
Smirnov (KS) statistics (IPress et al.ll200^ . The KS statis¬ 
tics looks at the maximum absolute difference between the 
empirical CDF of sample 1 and empirical CDF of sample 2. 
It is a nonparametric test that allows statistical comparison 
of two one-dimensional distributions. The KS test applied 
to the two samples yields 7)-value and P-value. Here, D is 
the maximum value of the absolute difference between the 
two CDFs. From Fig. [9] one can find that since the P-value 
is greater than our assumed significance level of 0.05, we 
reject the null hypothesis and conclude that the cumula¬ 
tive distribution of sample 1 is significantly different from 
sample 2 , with sample 1 being brighter and metal rich as 
compared to sample 2. Application of planefit procedure to 
the V-band data directly with distance determined using 
the MVF and MVC methods yields the value of c in Eqn.[9] 
to be —5.949 ± 0.233 kpc and —8.087 ± 0.364 kpc, respec¬ 
tively. Corresponding values of c obtained from the 7-band 
data calibrated to the V-band are —5.299 ± 0.509 kpc and 
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Table 1. Geometric parameters of the SMC determined fro m the OGLE III I- and V-band data using absolute magnitude values 
determined from the Fourier method and the relation given bv ICatelan &; CortesI l|2008h . 


Geometric parameters 


Band 


Method 

So/So 

Si/So 

S2/S0 


eionn 

MVF 

1.000 ±0.004 

1.334 ± 0.007 

08.884 ±0.123 

1.713 ±0.128 

74.063 ± 0.510 

MVC 

1.000 ± 0.004 

1.327 ± 0.008 

09.153 ±0.154 

1.912 ±0.136 

73.913 ±0.589 

MVF 

1.000 ± 0.002 

1.281 ±0.006 

07.067 ± 0.083 

2.860 ± 0.136 

72.848 ± 0.496 

MVC 

1.000 ±0.003 

1.243 ±0.008 

07.278 ±0.114 

3.637 ± 0.171 

72.516 ±0.714 











So/Mean(So) 


S,/Mean(So) 



Figure 6. Distribution of viewing angles and geometrical parameters of the SMC in histograms. Solid red lines denote Gaussian fit to 
the resulting distributions. 


—6.020 ±0.847 kpc. Negative values of c obtained in each of 
the four cases indicate that NE arm of the SMC is nearer to 
us than the plane of the SMC main body. Since in majority 
of the cases, we get the NE arm to be nearer to us by ~ 6 
kpc, we take this as the tentative estimate. It has also been 
found from the analyses that residuals of the planefit reveal 
a highly symmetric hyperbolic paraboloid warp of low am¬ 
plitude (~ 0.03 kpc). In order to test whether the results 
obtained so far are an artifact or a real f eature o f the S MC, 
we resort to the metallicity relation of ISmoled (l2005l l us¬ 
ing /-band data directly and absolute magnitude relation of 


ICatelan fc Cortej (l2008h as followed in section [9] Indeed, it 
has been found that same patterns are reflected and similar 
results are obtained in the resulting analyses. 


7 LINE OF SIGHT DEPTH OF THE SMC 

Line of sight (LOS) of the SMC has been studied us¬ 
ing numerous tracers in the literature and is an im- 
portant parameter t o ascertain the extent of the SMC 
jlCraM_et_^jJ_[200l|;_lSubramamm^&Subramania^ [200^, 

I 2 OI 2 I : Kapakos fc Hatzidimitrioijl2012l : Haschke et al.l[^12l l. 
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Table 2. Geometric parameters of the SMC determined fro m the OGLE III I- and V-band data using absolute magnitude values 
determined from the Fourier method and the relation given bv ICatelan &; CortesI l|2008h . 


Geometric parameters 




Method 

So/So 

Si/So 

S 2 /S 0 


01cnn 

SMC Main Body 

I 

MVF 

1.000 ±0.003 

1.185 ±0.004 

09.764 ±0.144 

0.753 ±0.093 

55.569 ±0.771 

(Ellipsoid-fit) 

MVC 

1.000 ±0.004 

1.186 ±0.005 

09.916 ±0.188 

0.573 ±0.106 

55.552 ± 0.949 


V 

MVF 

1.000 ±0.004 

1.346 ±0.007 

10.221 ±0.149 

0.801 ±0.106 

82.799 ± 0.401 


MVC 

1.000 ±0.006 

1.342 ±0.009 

10.839 ± 0.254 

0.354 ±0.106 

83.068 ±0.551 

SMC NE Arm 

I 

MVF 

— 

— 

— 

2.254 ± 0.856 

85.937 ± 14.030 

(Plane-fit) 

MVC 

- 

- 

- 

2.254 ± 0.864 

85.932 ± 14.149 


V 

MVF 

- 

- 

- 

1.284 ±0.210 

83.710 ±09.309 


MVC 

- 

- 

- 

0.500 ±0.193 

78.744 ±22.445 




Figure 8. Upper left panel shows separation of SMC into its northeast arm (shown as planar distribution) and that of the main body 
(shown as ellipsoidal distribution). The centre of the galaxy is marked by red star symbol. Other panels depict the projected distributions 
on the XY, YZ and XZ planes as well. The NE arm appears as a line in the two-dimensional YZ and XZ projections. 
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RR Lyrae stars in the present analysis have been used as 
proxies to get an estimate of the LOS depth. When combined 
with the metallicity values, it serves as a robust parameter to 
the underlying different galactic substr uctures llCrowl et al.l 
I 2 OOII : iKapakos fc HatzidimitriorJl2012l L To estimate the in¬ 
trinsic line of sight depth of the SMC, we take into account 
the standard deviation (cTobs) distribution of the RRab dis¬ 
tances and measurement errors ((Terror) in the individual dis¬ 
tances. The overall standard deviation (gobs) is a combin a- 
tion of these two factors in quadrature JCrowl et al. IH): 

^obs ■ '^int ^error• (10) 


The intrinsic ±lg line of sight depth is then given by 
(o'obs ~ 'terror) • Usiug the 7-band data, we have found the 
line of sight depth to be l(Jioa = 4.91 ± 0.65,4.12 ± 0.71 
kpc from the distances determined using MVF and MVC, 
respectively. Corresponding values of these parameters in 
the R-band were found to be 5.34 ± 0.61 and 4.63 ± 0.83 
kpc, res pectively. From a study of 454 R Rab stars in the 
R-band, iKapakos fc Flatzidimitri^ (l2012f l determined the 
LOS depth to be gmt = 5.3 ± 0.4 kpc. It may be noted 
that the LOS depth e stimate has a lower value when 
ICatelan &: Corti^ (120081 ) relati on is used to derive the line 
of sight distances then when dKovacs fc Walk^ 1200 ll ) re¬ 
lation is used. The variance weighted mean for the above 
estimates (weights = j-) y ields Oint = 4.81 ± 0.53 kpc 
dBevington fc Robinson! 2(303l) . There is a caveat in the LOS 
depth estimation from the optical band photometric light 
curve data as it is highly sensitive to the interstellar red¬ 
denings. A large sample of near infrared (NIR) photomet¬ 
ric light curve data can help us ascertaining the exact LOS 
depth of the SMC in different directions and thus help in 
locating the various components of the SMC. Use of differ¬ 
ent reddening maps applied to the optical band light curve 
data also yields differe nt estimates of the LOS depth val¬ 
ues l|Crowl et al.ll20f)ll ). However, most of the estimates of 
the Ig LO S depth values are limited in the range between 4 
and 6 kpc dCrow l et al.ll 2001 ^ Subramnman_&_Subramaniu^ 


l2012l : [Haschke et al.ll201^ Kanakos fc Hatzidimitriordl2012h . 

Despite completely different ways of obtaining the estimate 
of the LOS depth of the SMC from two OGLE bands, its 
value is consistent and well within the range. Combined with 
the LOS depth and the metalli city values of the sample, 
IKapakos fc Hatzidimitri^ (l2012l l found the existence of dif¬ 
ferent structures consisting of different populations. Divid¬ 
ing the sample on the basis of their metallicities, they found 
that the metal rich stars show small variations in their LOS 
depth while the metal poor stars exhibit larger LOS depth 
variations. From their analysis, they also concluded that the 
inner regions of the SMC consists of a thicker structure mim¬ 
icking a bulge. Metal rich and metal poor stars seem to 
belong to different dynamical structures. Metal poor stars 
have a smaller scale height and may belong to a thick disk 
and metal poor stars see m to belong to the halo. H o wever , 
the results obtained by IKapakos fc Hatzidimitri^ (l2012l 'l 
should be taken with caution since these have been obtained 
from a smaller sample of 454 RRab stars with their U-band 
light curves and ma y suffer from selection effec t . In o rder 
to test the results of IKapakos fc HatzidimitriorJ ||2012|) . we 
have also investigated the LOS depth of the sample using 
the 7-band data by dividing the sample into the metal poor 
([7^e/77] < [7^e/77] — g = —1.85 dex) and metal rich stars 


([7^e/77] > [7^e/77] -I- g = —1.60 dex). In order to find if 
there is any variation in the LOS depth of the two samples 
belonging to the different metallicity groups at a significance 
level of 0.05, we again perform a two sample KS test using 
the 7-band data calibrated to the U-band. The two sample 
KS test yields D = 0.110 and P-value = 0.90. The P-value 
> 0.05 indicates that we fail to reject the null hypothesis 
implying that there is no evidence in the data to suggest 
that the tw o CPFs are different. This i s cont rary to that 
obtained by IKapakos fc HatzidimitriorJ (120121 ). While cal¬ 
culating the LOS depth, it has been found that observed 
uncertainties in the distance measurement (gerr) of a few 
stars are greater than the observed depth (gobs)- The LOS 
depth determinations for those stars have been left out. 

To investigate the LOS depth variations in the SMC, 
we have divided the sample into three parts, viz., NE {x < 
—0.5, y > 0.5 kpc), central (—0.5 < x < 0.5,—0.5 < y < 
0.5 kpc) and SW (x > 0.5 ,y < —0.5 kpc). In order to see 
the variation of depth among the three different regions of 
the SMC, we study their normalized LOS CDF distribution 
through KS test. A two sample KS test indicates that that 
LOS depth of the NE and SW parts of the SMC are not 
significantly different. On the other hand, KS test applied 
on LOS CDF distribution of the central part and other parts 
combined (NE-I-SW) yields D = 0.316 and P-value=0.008. 
These values indicate that the difference in CDF of the two 
distributions are different at a significance level of 0.05, with 
LOS depth of central part being larger than rest of the parts 
of the SMC. This indicates that SMC may have a bulge. 


8 METALLICITY GRADIENT 

The presence of radial metallicity gradient in a galaxy in¬ 
dicates the presence of different stellar populations which 
are responsible for formation an d evolution of the galaxy 
(iBernard et al.l [20081 : ICionill2009l ). The detection of metal¬ 
licity gradient in the SMC usi ng the Call triplet (C aT) spec¬ 
troscopy was first reported bv ICarrera et al.l (l2008l ) using the 
350 red giant branch stars of 13 SMC fields distributed in the 
range 1° — 4° from the SMC center, wherein the metallicity 
decreases from the center of the galaxy towards the outer¬ 
most regions. In their search for the metallicity gradient they 
also found that this gradient is related to an age gradient 
where the stars concentrated in the central regions are gener¬ 
ally younger and metal rich. Recent studies attempt to bear 
out the presence of this gradient. Using the CaT metallicity 
study of the SMC field of 3037 r ed giant stars s pread across 
37.5 deg^ centred on the galaxy. [Pobbie et al.l (|2014|) found 
evidence of metallicity gradient of —0.075 ± 0.011 dex/deg 
over the inner 5 deg suggesting an increasing number of 
young stars with decreasing galacto-centric radius. However, 
the formal errors in their determinations of [7^e/77] are too 
optimistic. Finding meagre gradients as in the case of the 
SMC using the [7^e/77] values requires careful estimations 
of their errors and determining the weighted mean metal¬ 
licities and their errors while binning the galacto-centric ra¬ 
dius. The weighted mean metallicities and the errors in each 
bin tells how statistically significant the gradient is. Using 
a sample of 454 RRab stars in U-band from the OGLE-HI 
database, IKapakos fc Hatzidimitri^ (|2012| ) found a metal¬ 
licity gradient of —0.013 ± 0.007 dex/kpc with increasing 
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Figure 9. Cumulative distribution of extinction corrected RRab magnitudes Vq and metallicities [Fe/H] for the RRab stars (left panel 
shows the use of V-band directly on the empirical relation, whereas, right panel indicates the use of 7-band data calibrated to the 
V-band). NE arm and SMC main bosy are marked with red and blue color lines, respectively. 


metal abundance towards the dynam ical center of the SMC. 
However, the contrary was found by I Cion 1 ll2009ll using the 
[Fe/H] derived from the ratio of C- and M-type AGB stars 
analysed as a function of galacto-centric distance. A con¬ 
stant value of \Fe/H] = —1.25 ± 0.01 dex up to 12 kpc was 
obtained bv ICionil (l2009ll . The result of no metallicity gra¬ 
dient was reinforced by the study of the 25 SMC clusters 
on a homogeneous metallicity scale and with relative small 
metallicity errors bv IParisi et al.l (120091 '). Different values of 
metallcity gradients using various tracers obtained in the 
literature may be related t o their redistributions from th e 
positions they were formed llRoskar et al.|[2008l : ICion 1 1200911 . 
Comparison of metallcities of SMC held and cluster stars 
similarly may also yield signihcant metallcity gradients due 
to the fact that the held stars are more metal po or than ther 
corre sponding cluster stars surrounding them (IParisi et al.l 
l2009h . 

In the investigation of any gradient in a galaxy, it is 


important to determine viewing angles of the galaxy accu¬ 
rately. In this paper, we have determined the geometry of the 
SMC using the more precise value obtained from the princi¬ 
pal axis transformation applied on the de-projected Carte¬ 
sian coordinates. In order to investigate for any metallicity 
gradient in the SMC, the true galacto-cen tric distances were 
obtained using the procedure described bv ICion l (l2009ll . The 
coordinate system was rotated according to 

A'= Vcos(6i)-bTsin(6i) (11) 

V =-Asin(6l)-bycos(6(). (12) 

The {X', Y') system was then de-projected using 

Y” ^Y'/cos{i) (13) 

The angular distances were then calculated and converted 
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Figure 10. Metallicity distribution of RRab stars in the SMC. 
Metallicity values are binned on a 10 x 10 coordinate grid. In 
each bin, average metallicities and their associated errors are cal¬ 
culated which are due to uncertainties in the determination of 
Fourier parameter (pzi a^nd are shown in each grid box. 


Figure 11. Mean metallicity distribution of RRab stars as a func¬ 
tion of the galacto-centric distance (Rgc) with a bin size 

of 0.5 k pc. Mean metallicities in each bin have been obtained us¬ 
ing the lljurcsik & Kova^llOQOh empirical relations in the V and 
7-band with distances measured using the two methods described 
in the text. 


into kpc with 


ddeg = CX'^+Y"f (14) 

dkpc = D X tan (d), (15) 

where d is the angular distance to each star. Here the Carte¬ 
sian coordinates (X(q. d), y(Q, d)) are calculated fol lowing 
the method described in ivan der Marel fc Cionil ll200lll . This 
transformation helps in projecting a sphere onto a plane. 
Metallicity distribution of RRab stars in the SMC is shown 
in Fig. [TU] Metallicity values are binned on a 10 x 10 coor¬ 
dinate grid. In each bin, the average metallicities and the 
associated errors which are due to the uncertainties in the 
determination of the Fourier parameter (psi are calculated. 
Fig. [TT] depicts the variance weighted mean metallicity val¬ 
ues computed for four sets as a function of galacto-centric 
distance in bins of 0.5 kpc from the SMC center. A least- 
squares fit to the distance and mean metallicity values with 
their estimated errors to each of the data sets yields slopes of 
0.007 ± 0.038, 0.004 ± 0.040, 0.003 ± 0.064 and 0.006 ± 0.065 
dex kpc“^. All these values correspond to no statistically 
significant metallicity gradient. In the estimation of mean 
metallicity values in each distance bin, it has been ensured 
that the number of stars is greater than 10 for reliable statis¬ 
tics. Also, in the calculation of mean metallicity errors, the 
errors due to the uncertainties in the Fourier parameters 
and the systematic errors in each of the empirical relations 
have been taken into account. These two errors were added 
quadratically for each star in order to estimate the mean 
metallicity and its associated error in a distance bin. All 
of the empirical relations for metallicity calculations do not 
show any significant metallicity gradient with in the uncer ¬ 
tain ties, consistent with the results obtained bv ICion il (l2009l) 
and IParisi et al.l (|2009l) . 
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Figure 12. NE part of the SMC. Solid surface in green color 
denotes planefit to the data. 


9 SMC STRUCTURE DETERMINATION 
USING SMOLEC’S (2005) METALLICITY 
RELATION 

In order to compare, contrast and substantiate the viewing 
angle and geometrical structural parameter determinations 



referred to as ‘MVS’. The inclination angle and position 
angle of line of nodes for the entire sample of RRab stars 
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obtained from the ellipsoid fitting were found to be i = 
3M44±0°.129 and Sion = 74°.856±0°.514. The axes ratios 
were obtained as: 1.000±0.002 : 1.284±0.006 : 7.474±0.101. 
The sample of RRab stars are divided belonging to the SMC 
main body and NE arm on the basis of the changepoint 
analysis. Bimodal distribution of z-values also discriminates 
the NE arm from the main body of the SMC. A planefit 
procedure is applied to the NE arm which yields the value 
of c = —6.251 ± 0.854 kpc. The fitted plane is shown in 
Fig. 1121 The constant c represents the shift in the posi¬ 
tive or negative ^-direction from the xy-p\a.ne of the SMC. 
The value of c indicates that NE arm of the SMC is lo¬ 
cated at a distance ~ 6 kpc below the plane of the SMC 
main body, hence making it nearer to us by ~ 6 kpc 

than the main SMC plane. Several studies using Cepheids 
in the SMC found that the nor theast arm of the SMC is 
nearer to us than the central bar JCald well fc CoulsorJ[l98^ : 
iMathewson et al]|l988l : ISharpee et al.ll2002l l. This study us¬ 
ing statistically siginificant number of RRab stars further 
provides an insight into the issue thus confirming the pre¬ 
vious findings. The viewing angle parameters for the NE 
arm are i = 2°.205 ± 0°.407 and Sion = 85°.354 ± 6°.873. 
Fig. ini) shows the residuals after the NE arm planefit 
values were subtracted from their corresponding z-values. 
A highly symmetric hyperbolic paraboloid warp of very 
low amplitude 0.03 kpc) in the residuals is clearly dis¬ 
cernible. Tidal effects produced by the LMC (Large Mag¬ 
ellanic Cloud) and the SMC main body may be attributed 
to this warp. Similarly, the ellipsoid fitting procedure is ap¬ 
plied to the main body of the SMC which yields the values of 
i = 0°. 178±0°.091 and fton = 57°.053±0°.644 and axes ra¬ 
tios as 1.000±0.002 : 1.185±0.004 : 8.298±0.131. The pres¬ 
ence of metallicity gradient as a function of galacto-centric 
distance o f the SMC is a lso investigated using the empirical 
relation of ISmol^ (l2005l ) with galacto-centric radius deter¬ 
mined from inclination and position angle values as obtained 
from the J-band d a ta wh ere distance is determined from the 
ICatelan fc Corti3 (l2008l 'l empirical absolute magnitude re¬ 
lation. The gradient obtained is —0.008 ± 0.058 dex kpc“^. 
This is si milar to the gradient of 0.00 ± 0.06 dex kpc“^ ob¬ 
tained bv iHaschke et all (l2012l f corroborating evidence of no 
gradient and is con trary to the value of —0.013 ± 0.00 7 dex 
kpc“^ obtained by iKapakos fc HatzidimitriorJ (l2012l l from 
the analysis of OGLE-III E-band light curves of 454 RRab 
stars. Fewer number of data points and poor phase coverage 
in the case of E-band light curves of OGLE-III as compared 
to the 7-band light curves yield unreliable estimates of pa¬ 
rameters determined using them and hence should be taken 
with caution. The data in the E-band also suffer from spa¬ 
tial resolution which might also result into poor estimates of 
viewing angles and geometric parameter determinations for 
the SMG. Fig. [14] shows cumulative distribution of extinc¬ 
tion corrected RRab magnitudes Vo for the RRab stars in 
the calibrated E- band as well as metallicities calculated us¬ 
ing [Smole^ (l2005h relation. The probability values obtained 
from the two sample KS tests applied on stars belonging to 
the NE arm (sample 1) and the SMC main body (sample 
2) show that the CDFs of the two samples are significantly 
different at a level of 0.05 with sample 1 nearer to us and 
contains more metal rich stars as compared to sample 2, 
thus confirming the fact that the NE arm of the SMC is 
nearer to us and contains more metal rich stars. We have 



y [kpc] 

Figure 13. Distribution of residuals (z — 2^^) of the NE arm RR 
Lyrae stars as a function of y coordinates, where are the fitted 
2 -values obtained from Eqn. m . The distribution of stars take 
a hyperbolic paraboloid shape in 3d because of the gravitational 
interation exerted on it by the LMC and the main body of SMC. 
The scatter in the 2d plot is due to the projection from 3d. 


also investigated the LOS depth of the sample by divid¬ 
ing into the metal poor {[Fe/H] < [Fe/H] — a = —1.96 
dex) and metal rich stars {[Fe/F[] > [Fe/H] + a = —1.58 
dex). In order to find if there is any variation in the LOS 
depth of the two samples belonging to different metallic¬ 
ity groups at a significance level of 0.05, we again per¬ 
form a two sample KS test. The two sample KS test yields 
D = 0.150 and P-value = 0.79. The P-value > 0.05 in¬ 
dicates that we fail to reject the null hypothesis imply¬ 
ing that there is no evidence in the data to suggest that 
the two CDFs are d ifferent. This supports the result by 
iHaschke et al.l (l2012lf that there are no variations in the 
LOS depth for different metallicity groups along the different 
SMC fields as observed by OGLE III photometric survey and 
is con trary to that obtained by iKapakos fc Hatzidimitrioul 
||2012[) . While calculating the LOS depth, it has been found 
that the observed uncertainties in the distance measure¬ 
ments ((Terr) of a few stars are greater than the observed 
depth ((Tobs). The LOS depth determinations for those stars 
have been left out. The mean LOS depth (LOS) of the 
SMC is found to be 4.31 ± 0.34 kpc. This value is in good 
agreement with the values of 4.13 T 0.27 kpc, 4.0 7 ± 1.68 
kpc and 4.2 T 0.3 kpc found bv IKapakos et ^ ll201l]l. 
ISuf^manian fc SubramaniamI (I20l3 i and iHaschke et al.l 
(120121) , respectively. 

In order to see the line of sight depth variations in 
the SMC, we have divided the sample into three parts, 
viz., NE (a; < —0.5, y > 0.5 kpc), central (—0.5 < x < 
0.5, —0.5 < y < 0.5 kpc) and SW (x > 0.5 ,y < —0.5 kpc). 
In order to see the variation of depth among the three 
different regions of the SMC, we study their normalized 
LOS CDF distribution through KS test. A two sample KS 
test indicates that that LOS depths of the NE and SW 
parts of the SMC are not significantly different. On the 
other hand, KS test applied on the LOS CDF distribution 
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Vq [mag] [Fe/H](dex) 


Figure 14. Cumulati ve distributio n of extinction corrected RRab magnitudes Vq for th® RRab stars in the F-band with metallicities as 
determined from the llSmolec|l2005h relation. Stars in the northeastern arm and the SMC main body are marked with red and blue color 
lines, respectively. 


of the central part and other parts combined (NE+SW) 
yields D — 0.40 and P-valne=0.002. These values indicate 
that CDF of the two distributions are different at a sig¬ 
nificance level of 0.05, LOS depth of the central part be¬ 
ing larger than rest of parts of the SMC. Th i s sup ports 
the findings of ISubramanian fc SubramaniamI ll2009h and 
iKapakos fc Hatzidimitri~ I 2012h that the SMC may have 
a central bulge. 


Table 3. Parameters of exponential and King’s profiles for the 
radial number density distribution using the 4-band data 


Method 

n(roe) 

(stars kpc~^) 

n(rofc) 

(stars kpc~^) 

h 

(kpc) 

Tc 

(kpc) 

Tt 

(kpc) 

MVF 

844 ±73 

791 ± 62 

1.89 ±0.10 

0.72 ±0.10 

3.91 ±0.46 

MVC 

796 ± 69 

746 ± 60 

1.84 ±0.09 

0.72 ±0.10 

4.10 ±0.49 

MVS 

755 ± 67 

711 ±58 

1.81 ±0.09 

0.75 ±0.11 

4.11 ±0.49 


10 SCALE LENGTH AND TIDAL RADIUS OF 
THE SMC 


Scale length is defined as the point where the number density 
of stars falls by a factor of 1/e. Radial number density pro¬ 
files of RRab stars are obtained by finding the radial number 
density of stars around the centroid of the galaxy in a given 
galacto-centric radius with increasing binsize of 0.5 kpc and 
is shown in Fig llSI The number density profile is modeled 
with an exponential disk to find the scale length. To find 
the extent of the SMC in the xy-pVane, radial number den- 
sity distribu tion is fitted with King’s three parameter profile 
jKingll 19621) . The tidal radius gives an estimate of full extent 
of the SMC in t he xy-p\a.n e. King’s three parameter profile 
is described by llKindllo'^ ) 


n{r) = n(rofc) 


1 

[l + (r/r,)2]^ 


1 

[l + (n/r,)2]5 


2 


(16) 


The exponential disk profile is given by 

n(r) = n(roe)exp(-^), (17) 

where n(ro) and h represent the density of RRab stars near 
the galactic center and scale length, respectively, r is the 
galacto-centric distance in kpc. rt is called the tidal radius 
and represents the value of r at which n(r) falls to zero and 
Vc is the core radius where the number density of stars falls 


to half its central value JKin d ll962l) . The various parame¬ 
ters obtained from modeling the radial density profile with 
exponential disk profile and three parameter King’s profile 
are shown in Table [3] The variance weighted mean values 
were obtained as: = 0.73 ± 0.02 kpc, n = 4.03 ± 0.11 kpc 
and h = 1.84 ± 0.04 kpc. Although the values of h and 
found here are comparable to the values of 1.87 ± 0.10 kpc 
and 1.082 ± 0.02 kpc, the value of rt = 4.03 ± 0.11 kpc is 
drastically different from the v a lue of 7 ± 1 kpc obtained by 
ISubramanian fc SubramaniamI l|2012h . Tidal radius is com¬ 
parable to the la LOS depth of ~ 4.00 kpc obtained in this 
analysis. This implies that extent of SMC in the ly-plane 
is comparable to the front-to-back distance along the line of 
sight and points to the possibility that RRab stars in the 
SMC are distributed more like a spheroid/slightly ellipsoid. 
The values of rc and rt imply tha t the conce ntration pa¬ 
rameter defined as c = log (rt/rc) (iKin 3 ll962l) has a value 
- 1.71 ±0.07. 


11 SUMMARY AND CONCLUSIONS 

In this paper, we have studied both V and /-band light 
curves of RRab stars obtained from the OGLE-III project. 
Metallicities of RRab stars in the present sample were deter¬ 
mined using the l.Iurcsik fc Kovacd l|l99(j ) relation for both 
V and 7-band data. Absolute magnitudes of the RRab stars 
were determined in both ways from the Fourier parameter 
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Figure 15. Radial density distribution of the RRab stars in he 
SMC is shown as black circles. Error bars are due to the Pois- 
sonian distribution (cr y/N). Red and blue solid lines show 
the model fits with the empirical two param eter exponential and 
three parameter density profiles ijKin 3 ll962h . 


int er-relation of Kovacs fc Walk^ (1200 il l and those given 
by ICatelan fc CortesI (|2008l 'l. Results of the above analy¬ 


sis were furth er substantiated using the /-band data with 
ISmoled (l2005l ) metallicity rel ation for the /-band and empir¬ 
ical Mv — [Fe/H] relation of ICatelan fc CortesI (l2008l l. Also, 
since the /-band light curve data have very good phase cov¬ 
erage, we adopt the final results obtained from the /-band 
data only. Based on the metallicity of RRab stars combined 
with their distance distributions in the SMC, the following 
conclusions are drawn: 


(i) RR Lyrae stars in NE arm are nearer to us than the 
xy-plane of the SMC main body by a distance of ~ 6 kpc. 
Apart from that, the NE arm contains more metal rich stars 
as seen from metallicity distributions of these two groups 
using the two sample KS-statistics. Using the age, metal 
abundance and positio nal data of populo us cluster using 
Cepheids in the SMC, ICrowl et al.l (l200lh found that the 
eastern side of the SMC containing Cepheids is nearer to 
us and contains more metal rich stars but warrants further 
study. The investigation done in this study bears out this 
fact. The existence of an isol ated NE arm of the S MC was 
further corroborated using the lSchlegel et aP (Il998l l redden¬ 
ing map. However, it indicates that the NE arm is nearer to 
us by ~ 11 kpc than the galact ic plane of the SMC main 
body. The ISchlegel et al.l (Il998l ) maps use IRAS DIRBE 
data on the far-IR sky emission and estimate the extinction 
from the dust properties directly. Nonethless, the Schlegel 
maps are highly uncertain in the inner region of the SMC be¬ 
cause their temper ature structures were not sufficiently re - 
solved by DIRBE llSchlegel et al.l [19981 : IPessev et aDl2006lf . 
The choice of different search radius as well as different red¬ 
dening maps may have some effect on the final result of this 
paper. Availability of accurate and highly precise extinction 
measurements covering the entire region of the SMC will 


help enhance our understanding of its structure in great de¬ 
tails. 


(ii) Modeling the observed population of RRab stars in 
the SMC by a triaxial ellipsoid, we estimated the axes ra¬ 
tio, inclination of longest axis with the line of sight (i) and 
the position angle of line of nodes (0ion) on the sky from 
the analysis of the entire sample of RRab stars. Axes ra¬ 
tios of the galaxy were obtained as 1.00 ± 0.000 : 1.310 ± 
0.029 : 8.269 ± 0.934 with i = 2°.265 ± 0°.784, fton = 
74°.307 ± 0°.509 from the variance weighted /-band de- 
terminations . Usin g the 12 populous clusters in the SMC, 
ICrowl et al.l ll200lll obtained the axes ratios approximately 
as 1 : 2 : 4 modeling the SMC as a triaxial galaxy with 
declination, right ascension and line of sight depth as the 
three axes. On the other hand, using the same OGLE-II I 
SMC RRab dataset, ISubramanian fc SubramaniamI (l2012f) 
obtained the axes ratios, inclination of the longest axis with 
the line of sight (i) and the position angle of the pro¬ 
jection of the ellipsoid on the plane of the sky (0ion) as 
1.00 : 1.30 : 6.47, 0°.4, 74°.4, respe ctively, with no errors 
in the parameters quoted. Whereas, iHaschke et al.l (l2012|j 
found an inclination angle of 7° ± 15° and a position an¬ 
gle of 83° ± 21° from 1494 OGLE-III RR Lyrae stars using 
only the /-band data. These values are quite consistent with 
those obtained in the present study using distance distribu¬ 
tions determined from the Mv — [Fe/H] and Mv-Fourier 
parameters relations. 


(iii) NE arm of the SMC resembles a three-dimensional 
plane with the plane slightly warped in a hyperbolic 
paraboloid, while the the SMC main body is distributed 
like an ellipsoid. The NE arm has been htted by a plane 
equation to determine its viewing angles. The SMC main 
body has been htted with an ellipsoid obtained from the 
principal axis transformation of the moment of inertia ten¬ 
sor constructed from the {x, y, z) coordinates. Analyses in¬ 
dicate that the NE arm is slightly misaligned with the SMC 
main body, with a highly symmetric warp of low ampli¬ 
tude 0.03 kpc). The warp may have been caused by 
the tidal forces exerted on it by the main SMC body and 
the LMC. Axes ratios of the SMC main body were ob¬ 
tained as 1.00 ± 0.000 : 1.185 ± 0.001 : 9.411 ± 0.860 with 
i = 0°.507 ± 0°.287, 6»ion = 55°.966 ± 0°.814 from the vari¬ 
ance weighted /-band determinations. On the other hand, 
viewing angle parameters of the NE arm are found to be 
i = 2°.244 ± 0°.024 and 6»ion = 85°.541 ± 0°.332 from the 
variance weighted /-band determinations. 


(iv) Combining metallicities with spatial distribution of 
these tracers, no radial metallicity gradient in the SMC 
has been detected. Dividing the entire sample into three 
parts: northeastern (NE), central and southwestern (SW), 
we found that the central part has a signihcantly larger 
line of sight depth as compared to rest of the parts, indi¬ 
cating that the SMC may have a bulge. Independent anal¬ 
ysis done in this work exploring different empirical rela¬ 
tions and photo metric bands provides substantial evid ence 
to the findings of SuIam;Mniai 2 _fc SubramaniamI ll2009l ') and 
iKapakos fc Hatzidimitri~ I 2012h of SMC having a central 
bulge. 


(v) Radial number density of RRab stars in the /-band 
were model ed by an exponential disk and three parameter 
I Kind lll962ll profiles. Following values were obtained from 
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the modeling: h = 1.84 ± 0.04 kpc, = 0.73 ± 0.02 kpc, 
rt = 4.03 ± 0.11 kpc and c = 1.71 ± 0.07. 
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